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Abstract 

There is a good deal of current interest in the condensate < g^A'^^Aa > which has recently be 
shown to be the Landau gauge version of a more general gauge-invariant expression. In the present 
work we consider quark propagation in the presence of such a condensate which we assume to 
be present in the vacuum. We describe the vacuum as a random medium of gluon fields. We 
discuss quark propagation in that medium and show that the quark propagator has no on-mass- 
shell pole indicating that a quark cannot propagate over extended distances. That is, the quark is 
a nonpropagating mode in the gluon condensate. 
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I. INTRODUCTION 

Some years ago we presented a discussion of the properties of the gluon condensate [1,2]. 
Central to our work was the specification of the condensate < g'^A'^ > which was obtained 
from the known value of < | {(Xs/t^)G'^^G'^^^ | >. (The value of the latter quantity was 
the only parameter we needed for our analysis.) We obtained a dynamical gluon mass of 
mG=649 MeV and a dynamical quark mass of m^'=432 McV, which arose from coupling to 
the gluon condensate [1,2]. In addition, we predicted a glueball mass of 1.0 ± 0.15 GeV. 

Our work was limited by the fact that the condensate < g^A^ > was thought to be non 
gauge invariant. However, recent years have shown a renewed interest in such a condensate. 
For example, we note that the operator G'^^^G'^^ generates terms of order in the operator 
product expansion, while Afj^A^ is the only term that has a vacuum expectation value which 
will generate l/p^ terms in studies making use of the operator product expansion. There is 
significant empirical evidence for the importance of such terms [3]. 

The issue of gauge invariance has been discussed by several authors [4,5] and it was argued 
that < > may be important for the study of the topological structure of the vacuum 
and quark confinement [4,6,7]. Kondo [4] was responsible for introducing a BRST-invariant 
condensate of dimension 2 



where c{x) are Faddeev-Popov ghosts, a is the gauge-fixing parameter and Q is the inte- 
gration volume. To quote Kondo: "It is clear that O reduces to in the Landau gauge 
a = 0. Therefore, the Landau gauge turns out to be a rather special case in which we do 
not need to consider the ghost condensation." Here, the minimum value of the integrated 
squared potential is which has a definite physical meaning [4]. 

Recently Arriola, Bowman and Broniowski [8] have made use of lattice data for the quark 
propagator. They have considered an operator product expansion (OPE) in the analysis of 
the QCD lattice data, including < A^ > jQ"^ terms and use the relation mQ — {?,/'i2)g'^ < 
A^ > to obtain mc = (625 ± 33) MeV, a numerical result that is consistent with the one we 
have obtained by a different method [1,2]. 

Our work contains the following material. In Section II we will introduce a description 
of the QCD vacuum as a random medium and define our model for the matrix elements of 




(1.1) 
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the condensate field. In Section III we review the theory of wave propagation in a random 
medium and introduce the quark propagator. In Section IV we calculate the quark self- 
energy, E(p) = A{p'^)'^ + BIp^), and present results for M{p^) = + mg"'')/(l — 
B{p^) + nf^"^ and [1 — in Figs. 2-14. Section V contains some further discussion and 
conclusions. 



II. THE QCD VACUUM AS A RANDOM MEDIUM 

We now proceed to consider quark propagation in a random medium which is character- 
ized, in part, by finite value of the < > condensate. It is useful to separate the vector 
potential into a condensate field and a fluctuating field. We write 

A^(x)=A«(x)+A"(^)- (2.1) 

This division of the vector potential into a condensate field, which we assume may be treated 
as a classical random field, and a fluctuation field, A'^^{x), will be useful. The basic idea 
is that, independent of the nature of the condensate, be it a monopole vortex or instanton 

condensate, there is a natural separation of the field into a nonperturbative condensate 
field and a fiuctuation field which can be treated perturbatively. (Some discussion of this 
separation is given in Shuryak's book with reference to its use in a QCD sum-rule analysis 
[9].) 

We introduce a vacuum state which has the following properties: 

< vac I A'„(0) I vac >= 0, (2.2) 

< vac I AUO)A^,(0) I vac >= ^f<i>ff, (2.3) 

etc. 

We will find it useful to generalize this result to provide a covariant description of the 
vacuum state. We define a vacuum state, | vac >, which has the following properties 

< vac I A^(0) \vac>=0, (2.4) 
< m I /A«(0)At(0) I ^c >= -^-fg'cl^ff, (2.5) 
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< vac I A;^(0)A^(0)A^(0) I mc >= 0, (2.6) 

< vac I /A;^(0)A^(0)A^(0)AJ(0) I vac > (2.7) 
^<vac\ A^(0)A^(0) \vac><vac\ A^(0)A^(0) | vac > 
+ <vac\ A^(0)A^(0) \vac><vac\ A^(0)Aj(0) | vac > 
+ <vac\ A)^(0)A^(0) \vac><vac\ A|^(0)A^(0) | vac > 

= {gixi^gpri^ab^cd + gtipgvn^achd + gmgi^pSadhc)g^(Pt/ (16 • 64), (2.8) 

etc. Here, all odd correlation functions are taken to be zero and all even correlation functions 
may be expressed in terms of the two-point correlation function, as in Eq. (2.8). This 
characterization is that of a Gaussian random medium. 
More generally, we might put 

< ^c I g'A;{x)Al{y) \ ^c >= -^g'^f-^exp [- \ {x - yf /L] (2.9) 

where L is a correlation length. The model we will study has L — > oo. That is, L will be 
taken to be large compared to the characteristic mean-free-path associated with the damping 
of the quark propagator. 

We wish to study the equation for the quark field 

{i^^d, - gK{x)i,^)ij{x) = (2.10) 

and the associated quark propagator. We have only included the condensate field in Eq. 
(2.10) and will treat A^(a;) as a random field having the properties described above. With 
that goal in mind, we will review the classical theory of wave propagation in random media. 

III. WAVE PROPAGATION IN A RANDOM MEDIUM 

A useful review of wave propagation in random media has been given by Prisch [10]. We 
follow the notation of that work, with a few minor exceptions. One of the more elementary 
problems, which is characteristic of the statistical analysis, is that of scalar wave propagation 
in a time-independent, lossless, homogeneous, isotropic random medium. The field satisfies 
the Helmholtz equation 

VV(r) + kln^{r)^{f) = (5(r), (3.1) 
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where n^(r) is the index of refraction and ko is the free-space wave number. One writes [10] 

n'(r) = l + //(r), (3-2) 

where /x(r) is a centered homogeneous and isotropic random function with finite moments. 
That is < //(r) >= 0, while the correlation function 

r(| f - P I) =< n{r)n{r') > (3.3) 

is unequal to zero. Here the brackets denote an ensemble average. For simplicity, we can 
assume that ^{r) is a Gaussian random function. 

The Green's function satisfies a random integral equation. With 



47r I f — r' I 

we have 

G{r,r') = G(°)(f,r') - J G^'^\f,r-{)ii{fl)G{fl,?)dfl. 



(3.5) 



We write the last equation as 

G = G^°^ - G^^^LiG, (3.6) 

where Li is the random element. We are interested in the ensemble average of G, which we 
denote as < G >. Thus 

< G >= G^^)- < G^^^LiG >, (3.7) 

since < G^'^^ >= G^^\ (Here /cq has been absorbed in Li.) 

One may develop a diagrammatic analysis of a perturbative expansion of Eq. (3.7) [10]. 
When averaging the various terms in such an expansion over the ensemble one finds it useful 
to introduce so-called p-point correlation functions rp(l, 2, ••• ): 

ri(i) = <//(i)>, (3.8) 

= 0, (3.9) 
r2(r,2)=<//(l)//(2) >, (3.10) 
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r3(l,2,3)=<//(l)//(2V(3)>, (3.11) 
r4(r,2,3,4) =< /x(r)/i(2)/i(3)/i(4) > -12(1,2)12(3,4) (3.12) 

-r2(r,4)r2(2,3)-r2(r,3)r2(2,4), 

etc. For a Gaussian centered random function Fi = 0, F3 = 0, F5 = 0, etc. 
A simple first-order smootliing approximation is 

< G{r, ?) >= G(°)(f, ?) + k^J G(°)(f, 1)G(°)(1, 2)F(r, 2) < G(2, ?) > dld2, (3.13) 

in the notation of Ref. [10]. (See Eq.(4.33) of [10].) This approximation is most appropriate 
for a small (generalized) Reynolds number defined as 

R = eklL, (3.14) 

where e is some (dimensionless) measure of the scale of /i [/x = e/i*, where /i* = 0(1)]. 

We will here be interested in the case of large generalized Reynolds number. In this case 
Eq. (3.13) is replaced by the Kraichnan equation [11] 

< G{f, ? >= G'(°)(f, P) + k^J G'(°)(f, 1) < G(l, 2 >< G(2, P) > F(l, 2)dld2, (3.15) 
or in a more abstract notation 

< G G^") + < Li < G > Li >< G > . (3.16) 

(See Eqs.(4.76) and (4.79) of Ref. [10].) Kraichnan [11] has shown that Eq. (3.7) may be 
reduced to Eq. (3.16) upon the introduction of an additional stochastic element (random 
couphng between different Fourier components of the field). It has also been shown [11] 
that Eq. (3.16) will give sensible results, since it can be considered both as the solution of 
a model equation and as approximate solution of Eq. (3.7). 

It is useful to introduce a self-energy operator, E. From Eq. (3.15), we have 

< G(f, P) >= G(°) (f, P) +kl j G(°) (f, 1)S(1, 2) < G(2, P) > dld% (3.17) 

where 



E(l,2) =< G(1,2)>F(1,2). 
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(3.18) 



FIG. 1: Equations for the ensemble- averaged quark propagator and the self-energy in the random- 
coupling model, a) The solid line is the full propagator and the light line is the propagator 
[p — m^^]^^. b) A self-consistent model for the (irreducible) self-energy 

S(p2). The cross-hatched 

circle denotes the condensate and the wavy lines are condensate gluons of zero momentum, c) Here 
we show the type of nested diagrams which are summed in this model for the self-energy S(p^). 

As usual, this equation is particularly simple in momentum space, where we can write 

" ""^"^ |G(°)(rti---s(p) - '"■I"' 

in the case of an infinite medium. The equations for the propagator and for S(p) are shown 
in Fig. 1. There we also show the kind of diagrams which are being summed in the random- 
coupling approximation. 

The translation of this analysis to a study of the quark propagator is immediate. We 
define a free propagator, S^^\x, x'), and the ensemble average of the correlation function of 
the Heisenberg fields: 

< S{x,x') >=< vac I T(-?/'(a;)^/'(x')) | vac > . (3.20) 



The random element, Li, is 

Li = yA^(x)7;.A72. (3.21) 

We write 

rj(l, 2) = < ^ I /A^(0)A5(0) I > (3.22) 

Here we have generahzed the scalar correlation function, r(l,2), to a form appropriate in 
the case the correlation functions involves random functions with color and Lorentz indices. 
In momentum space we have 

S^'\p)^l, (3.24) 
P 



and 



<SiP)> = J^y (3-25) 

We need not include an ie in the denominators on the right-hand side of Eqs. (3.24) and 
(3.25) since we will consider solutions for < S{p) > which have no singularities for real p^. 
This feature of the analytic structure of the propagator is analogous to that found when 
studying classical wave propagation in a random medium [10], as noted earlier. 

IV. THE QUARK SELF-ENERGY IN THE KRAICHNAN RANDOM- 
COUPLING APPROXIMATION 

The approximation considered here is shown in Fig. 1. There the cross-hatched region 
denotes the condensate and the wavy lines are condensate gluons. When translated into 

momentum space, the assumption of large correlation length (L — >■ oo) means that the 
condensate gluons may be taken to carry zero momentum. 
We write the quark self-energy as 

S(p) = A(p2)|5 + 5(/). (4.1) 
In the random-coupling approximation we have [See Fig. 1], 



where m^"'' is a "current" quark mass and = g^(t)Q/2A contains a color factor of 4/3. 

Let us first consider the case m^"^ = 0. We see that Eq. (4.2) exhibits a bifurcation 
phenomenon. There is a solution with S = 0; however, we will consider another solution 
with i? 7^ 0, which breaks chiral symmetry. That solution is 

= (4.3) 

B{p^) = 2\J'p^ + forp^ > -h^. 



and 



A{p^) = ^ [l - Vl - 8«7p'] , (4-4) 



S(p2) = 0, for p2 < -K^. 
It is useful to define a dynamical mass 



which for rrf^^ = is 



M(p2) = (p2 + «;2) 1/2^(^2 ^ ^2)_ (4_g) 

It may be seen that Eqs. (4.3)- (4.4) represent the m™'' = limit of a continuous solution 
of Eq. (4.2) obtained with m^"'' ^ 0. Equation (4.2) may be solved for m™'' ^ 0. One finds 
that Mij?) satisfies a fourth-order polynomial equation. That equation has four solutions 
for each p^. We choose a real continuous solution such that M(p^) ■nf^'^ as — oo. 

The value of M^(j9^), B{p^) and [1 — A(p^)] for such solutions are shown in Figs. 2 to 14, 
for various values of m^'"'. The values chosen are appropriate for up, down, strange, charm 
and bottom quarks. [We have also taken k^=(232 Mev)^.] 

We note that for the solutions chosen, for large spacehke we have A{p'^) — > 0, S(p^) — > 
and 

< Sip) >^ j-^,. (4.7) 
P — m^^^ 

For timelike p"^ we have MijP') p^ + for large p^. It may be seen that the equation 
p^ — M'^{p^) has no solution and therefore there are no poles in the statistically averaged 
quark propagator for realp^. This generalizes the corresponding result of the classical theory 
to the relativistic theory. 
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FIG. 2: The square of the dynamical mass, M^(p^), for up and down quarks. Here we chose 



m„ 



=5 MeV. [For > 0, M^{p^) = + k^.] 




FIG. 3: The mass parameter [B{p'^) + m^"*"] for up and down quarks with m^"'' =5 MeV. 
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FIG. 4: The wave function renormalization parameter [1 — A{p'^)] for up and down quarks. 




FIG. 5: The square of the dynamical mass, M^lp"^), for the strange quark. Here 171^^^=125 MeV. 
(Note the change in scale with respect to Figs. 2-4.) 
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FIG. 6: The mass parameter [-B(p^) + m™''] for the strange quark. (Here m^"'' = 125 MeV.) 
V. DISCUSSION 

In this work we have shown that the ensemble average of the quark propagator exhibits 
damping. A similar result was obtained for the gluon propagator in an earlier work [12] 
where we used the first-order smoothing approximation [13] to generate an expression for the 
vacuum polarization operator. In Ref. [12] we studied the gluon propagator in momentum 
space and carried out a Fourier transformation to coordinate space. In that case one could see 
explicitly how the absence of on-mass-shell singularities in the propagator leads to damping 
of the wave both for spacelike and timelike propagation. In that analysis, carried in the 
Landau gauge, it was seen that the gluon obtained a dynamical mass uiq = g'^(t)l/4:. We 
also note that in a lattice simulation of SUc{5) Yang-Mills theory in the Landau gauge [14] 
it was found the the gluon obtained a dynamical mass of mo = 600 ± 90 MeV. Since we had 
put = we use this information to fix the value of k [k—232 MeV] . 

In the work reported here we find that the solution which exhibits nonpropagation of 
quarks over large distances is also characterized as having B 0. Nonzero values for B 
denote chiral symmetry breaking. It is seen that only a single parameter, g'^cj)'^, governs 
both phenomena. This is a satisfactory result in the case of QCD, since only a single 
dynamically generated dimensionful parameter should characterize the theory. 
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for the strange quark. [See 
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FIG. 7: The wave function renormalization parameter [1 — 
caption to Fig. 5.1 
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FIG. 8: The function A^p^) for the strange quark (m^"^=125 MeV) for an extended range of values 
of p"^. 
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FIG. 10: The mass parameter [B{p'^) +171^^^] for the charm quark. (Here m^"*" = 1375 MeV.) [See 
caption to Fig. 9.] 



nr>(f"' = 1375 MeV 




FIG. 11: The wave function renormahzation parameter, [1 — A[p'^)] for the charm quark. Here 
^cur _ ]^375 GeV. (Note the change of scale with respect to Figs. 2-8.) 
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FIG. 12: The square of the dynamical mass, M^{p'^), for the bottom quark. Here m^"'' = 4750 
MeV. (Note the change of scale with respect to Figs. 2-11.) 



15 




__| 1 i ,| ,r ti- J.„„A„,fli-.^ ^ # 


i.. 



FIG. 13: The mass parameter [B(p^) + m'^^^] for the bottom quark. Here m'^'^ = 4750 MeV. [See 
caption to Fig. 12.] 
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FIG. 14: The wave function renormalization parameter [1 — ^(p^)] for the bottom quark. Here 
m'q'''' = 4750 MeV. [See caption to Fig. 12.] 
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